The Lower Critical Dimension of the XY Spin Glass 
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We investigate the XY spin-glass model in two and three dimensions using the domain-wall 
renormalization-group method. The results for systems of linear sizes up to L = 12 (2D) and 
L — 8 (3D) strongly suggest that the lower critical dimension for spin-glass ordering may be d c ~ 3 
rather than four as is commonly believed. Our 3D data favor the scenario of a low but finite spin- 
glass ordering temperature below the chiral transition but they are also compatible with the system 
being at or slightly below its lower critical dimension. 



PACS numbers: 64.60. Cn, 75.10.Nr, 05.70.Jk 

It has been known since the early work of Villain 
that frustrated planar models possess, in addition to the 
continuous degeneracy associated with spin rotations, a 
discrete two-fold degeneracy associated with the invari- 
ance of the Hamiltonian under a reflection about an ar- 
bitrary direction. As a consequence, each plaquette has 
an Ising-like degree of freedom, the chirality, that cor- 
responds to the handedness of the configuration of the 
spins around it. It is also well established that in 
two and three dimensions chiral and spin variables de- 
couple at long distances and order separately Q. In 
2D long range chiral and spin order appear simultane- 
ously at zero temperature, but the transition is unusual 
in that there are two independent divergent correlation 
lengths as T -► |,|. While there is convincing evi- 
dence that in 3D the chiral-glass transition occurs 
at a finite temperature Tcg, existing Monte Carlo simu- 
lations |5|,|7j and studies of the scaling of defect energies 
at zero temperature suggest that spin-glass order 
sets in at Tsg = just as in two dimensions. These, 
as well as older numerical results |§|-|i"c|] have led to the 
belief that the lower critical dimension (LCD) for spin- 
glass order in this model is d c > 4, a conjecture that is 
widely accepted even if a rigorous proof has turned out 
to be elusive [pd|-[l^1 . The limitations inherent to the nu- 
merical methods raise some doubts about the robustness 
of this conclusion, however. The systems studied with 
the defect-energy method are rather small. While the 
domain wall energy flows towards weak coupling with in- 
creasing linear size L for systems with L < 4 p|, more 
recent simulations pj for L < 6 show a tendency towards 
saturation for the bigger sizes. It is thus conceivable that 
the systems studied up to now are still far from the scal- 
ing regime. The Monte Carlo simulations were for 
much larger lattices but the rapid increase in the ther- 
malization times with decreasing T makes it impossible 
to attain the relevant temperature region, T <C Tcg- 
It appears that is also difficult to reach the asymptotic 
regime with this technique. Indeed, Kawamura p] has 
noticed that the Binder function associated with the spin- 



glass order parameter fails to show scaling behavior in the 
temperature range covered by his simulations, T > Tcg- 
The spin-glass susceptibility does seem to scale but Jain 
and Young [Q found that their data can be fitted with 
comparable accuracy by assuming such different values 
of the spin-glass critical temperature as Tsg = and 
Tsg = 0.45 J. In view of such uncertainties we may still 
regard the nature of the low-temperature phase of the 3D 
XY spin-glass as an open problem. 

In this paper we study the XY spin-glass model 
in two and three dimensions with the domain-wall 
renormalization-group method (DWRG) §§ . Using a 
new and powerful algorithm for the search of ground- 
state energies we have been able to study systems with 
linear sizes up to L = 12 (in 2D) and L = 8 (in 3D). In 
both cases our largest system contains more than twice 
as many spins as have been considered before. Our re- 
sults in 2D agree with those found by other authors [0, 
confirming that the scaling regime had been reached in 
previous simulations. This turns out not to be the case in 
three dimensions, where we find that there is a crossover 
between small- and large- T behavior at L rj 5. The 
domain- wall energies scale as W ~ L x . The sign of the 
stiffness exponent A may be positive or negative depend- 
ing on whether the system is above or below its LCD. We 
find A c = 0.5 ± 0.2 and A s = 0.056 ± 0.11 for chiral and 
spin domain walls, respectively. The fact that A c > 
is the signature that there is long-range chiral order at 
low-temperature in the system as found by other authors 
P,||. The smallness of |A S | constitutes to our knowledge 
the first numerical evidence that the LCD of the (±J) 
XY model may be close to three. Two scenarios are 
compatible with the measured value of A s , i) a spin-glass 
transition at a finite temperature Tsg *C Tcg or ii) a 
zero-tempcraturc transition with the correlation length 
£ ~ T~ v and v ^S> 1. Statistics favors the former. 

The Hamiltonian of the model is 



H = -J2 J v cos (Oi-Oj), 



(1) 



where the sum runs over all pairs of nearest-neighbor 
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sites of the 2D square or 3D cubic lattices. The ex- 
change couplings Jij are random independent variables 
that take the values J = +1 and J = — 1 with equal 
probability. In the DWRG method ||,|)| one studies the 
sensitivity of the system to changes in the boundary con- 
ditions at T = 0. The ground-state energies of an ensem- 
ble of systems of size L D are calculated using periodic 
(P) and anti-periodic (AP) boundary conditions along a 
given direction while keeping fixed boundary conditions 
along the [d — 1) remaining directions. The width of 
the distribution of differences of ground-states energies, 
W S (L) = [(Ep -Ea) 2 ] 1 / 2 is interpreted as an effective 



coupling constant between blocks of L spins. This is 
expected to scale as W(L) ~ L x for large enough L. 
If the stiffness exponent, A, is positive the rigidity of a 
block diverges with its size, a sign that there is long- 
range order in the system. If A < there exists a length 
scale L~I A I ~ T beyond which the effective coupling 
between blocks becomes smaller than the temperature. 
This scale is identified with the correlation length and the 
correlation-length exponent is obtained from v = 1/|A 
p0| . Chiral ordering may be studied similarly except 
that the calculation requires knowledge of the ground- 
state energy for reflective boundary conditions as well as 
the other two 

The success of this approach relies on the availability 
of an efficient and accurate way of finding the lowest- 
lying states of the system. In the usual spin-quench 
algorithm 14 (SQA) one randomly generates long se- 
quences of metastable configurations among which one 
hopes to find the ground-state or states sufficiently close 
in energy. The number of metastable states of a frus- 
trated system increases rapidly with its size (probably 
exponentially jnj ) and so does the number of trials that 
need to be performed to have a non-negligible chance 
of finding a relevant state. This fact restricts severely 
the sizes of the systems that can be studied using this 
method. In order to go beyond the limits of the SQA 
we have developed a ground-state search algorithm of far 
greater efficiency. This algorithm is inspired by some of 
the morphological features of the low-lying states of the 
frustrated XY model that are as follows jl5|,[l6]]. i) In a 
given sample there are regions with a low density of frus- 
trated plaquettes where the spins form almost collincar 
domains, and regions where frustration is high and the 
spin arrangement looks random. The position, size and 
shape of the domains are mostly determined by the ran- 
dom bond realization and hardly vary from one state to 
another, ii) Apart from smooth spin-wave-like distor- 
tions, the essential differences between the spin configu- 
rations of any two low-energy states are large amplitude, 
almost rigid rotations of the individual domains, and chi- 
rality reversals of a fraction of the frustrated plaquettes 
in the regions separating them. The energies of states 
that do not have this structure are much higher p6|,p^[. 



Our method suceeds in preserving the domain structure 
at each stage of the procedure by treating the spins in 
high local fields (domain spins) differently from those in 
the more frustrated regions. In this way the appearance 
of uninteresting high-energy configurations becomes un- 
likely. The algorithm will be sketched in the following 
paragraph. Full details will be given in a separate publi- 
cation Jl6[ . 

The initial (or parent) state of the sequence, {0q}, is 
obtained by a conjugate-gradient minimization (CGM) 
of the energy with a random spin distribution as ini- 
tial condition. Next, new spin configurations are gener- 
ated by iterating the following loop any number of times, 
i) The spins in the n-th configuration {9 n } are divided 
into two groups according to whether their local field 
hi = Y]j cos (9i — 9j) is greater or smaller than an 
appropriately chosen threshold field, < h t < 2D (see 
below). The spins in the first group are defined as the 
domain spins, ii) Existing correlations between the do- 
mains and the rest of the system are broken by means of 
a random rigid rotation of the former, iii) A fraction p 
of the spins in weak local fields are randomly picked and 
their orientations are randomly reset, iv) The energy of 
the subsystem of domain spins is minimized with the rest 
of the spins held fixed, v) The resulting configuration is 
let to relax performing a CGM of the total energy of the 
system. The outcome is the next state in the sequence, 
{#n+i}. vi) The energy of this state is stored and h t is 
rescaled if necessary (see below), vii) End of the loop. 
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FIG. 1. Distribution of the energy per spin E of two series 
of 5 000 metastable states each of one realization of the XY 
J — ±1 spin glass on a 512-site 3D lattice. The states have 
been generated using the spin-quench method (a) and the 
algorithm described in the text (b). 

The parameter h t must be carefully chosen for this al- 
gorithm to be efficient. It defines the degree of collinear- 
ity required of an ensemble of spins for it to qualify as a 
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domain (h t = 2D for perfectly aligned spins in an homo- 
geneous environment). If h t is set very high, too many 
spins are involved in step iii) above and the domain struc- 
ture is not preserved. In this regime the method is equiv- 
alent to the SQA where all the spins are randomly reset 
at each step. On the other hand, too small a thresh- 
old field leads to trapping, i.e., all the states of the se- 
quence are in the vicinity of the parent state. There is 
an optimum value of the threshold field in between these 
extreme cases. In the running version of our code the per- 
formance of the algorithm is continuously monitored and 
a procedure has been devised that allows h t to self-adjust 
when degradation is detected. This method is particu- 
larly appropriate for the treatment of large systems. This 
is illustrated in Fig. [I] which shows the energy distribu- 
tions of two sequences of metastable states of the same 
realization of a disordered 3D system of N = 8 3 spins 
with anti-periodic boundary conditions. Histograms (a) 
and (b) have been obtained from 5000 energy levels each 
generated using the SQA and our algorithm, respectively. 
In the latter case we have verified that trapping had not 
occurred by making sure that the same histogram results 
from sequences originating from different parent states. 
Distributions (a) and (b) are quite different. Whereas the 
histogram of the energy levels obtained with the standard 
method is wide and peaks at high energies, that of our 
sequences is much narrower and is concentrated in the 
low-energy end of the spectrum. It is remarkable that 
most of the energies found with the new method lie in 
a region where the SQA did not find any state after the 
same number of trials. The calculation with our method 
takes only 20% more CPU time than that with the tra- 
ditional one. 

We have used this method to perform DWRG calcu- 
lations for two- and three-dimensional systems. Much is 
known about the two-dimensional case which therefore 
serves as a test of our methods. We have determined 
spin and chiral defect energies for systems of size £ x £ 
with £=4,5,6,7,8,10 and 12. The disorder averages were 
taken over 128000 (£=4), 64000 (£=5), 12800 (£=6 and 
7), 6400 (£=8), 2560 (£=10) and 1280 (£=12) indepen- 
dent bond configurations, respectively. The ground state 
energy was estimated from the analysis of sequences con- 
taining 20 (£=4), 50 (£=5), 100 (£=6), 200 (£=7 and 8), 
1500 (£=10), and 3000 (£=12) states. Low-energy states 
were accepted as ground-state candidates only if they and 
their chirality-reversed partners appeared several times 
in the sequence at widely spaced positions. This guar- 
antees that the states in the sequence come from well 
separated regions of phase-space. A log-log plot of the 
size dependence of the two-dimensional defect energies is 
shown in Fig. |^. The symbols are the numerical results 
and the dashed lines are fits to a power-law. The fits 
are of good quality even for small £. The renormalized 
stiffness decreases with increasing £ for both types of do- 
main wall indicating that spin and chiral variables only 



order at zero temperature. From the slopes determined 
by the fits we can compute the correlation length expo- 
nents for the two transitions. We find v s = 1.29 ± 0.02 
and v c = 2. 57 ±0.003 for the spin and chiral order param- 
eters, respectively. These values are in good agreement 
with the results of Kawamura and Tanemura Q who find 
v s = 1.2 ± 0.15 and v c = 2.6 ± 0.3, respectively in their 
DWRG calculations for systems with £ < 8. Moreover, 
the value of the chiral exponent deduced from our data 
is very close to the correlation length exponent of the 2D 
Ising spin-glass obtained by Monte Carlo |l7]] (2.6 ± 0.4) 
or transfer-matrix [|l8| (2.59±0.13) methods. This agrees 
with the view that the chiral transition in the XY spin 
glass model is in the Ising universality class. 
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FIG. 2. The L-dependence of the spin (squares) and chi- 
ral (circles) domain-wall energies of the (±J) XY spin-glass 
model on L x L square lattices. Where not drawn the statis- 
tical error bars are smaller than the size of the data points. 
The dashed lines are the power-law fits mentioned in the text. 

We now turn to the results for three-dimensional sys- 
tems. Much longer sequences of stationary states are 
needed to determine the ground-state energy in this case. 
Our series consist of 100 (£=3), 200 (£=4), 500 (£=5), 
1000 (£=6), 3000 (£=7) and 5000 (£=8) states, re- 
spectively. Moreover, the results for the bigger lattices 
were further checked by repeating the calculation start- 
ing form three different parent states. Sample averages 
were taken over 25600 (£=3), 6400 (£=4), 1280 (£=5,6), 
640 (£=7) and 128 (£=8) samples, respectively. The 
£-dependence of the effective couplings in three dimen- 
sions is shown by the symbols in Fig. |^. An important 
difference with the 2D case is that a crossover region is 
clearly seen around £=5. Power-law behavior, shown in 
the figure by the dashed lines, is only observed for the 
largest lattices. We first discuss the chiral degrees of 
freedom. As found in previous studies @||, W c flows 
towards strong coupling signaling the existence of finitc- 
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temperature chiral-glass transition in the infinite system. 
From the fits shown in the figure we can estimate the chi- 
ral stiffness exponent A c = 0.56 ± 0.18. 




FIG. 3. The L-dependence of the spin (squares) and chi- 
ral (circles) domain-wall energies for the (±J) XY spin-glass 
model on L x L x L simple cubic lattices. The dashed lines 
are the power-law fits mentioned in the text. 

The spin domain-wall-energy decreases rapidly with L 
for small sizes but it exhibits a much slower variation for 
L > 5. This behavior was not observed in older simu- 
lations Q but Kawamura's more recent DWRG results 
distinctively show the beginning of the saturation of the 
spin defect energy for L ~ 5. Fitting the results for 
the four largest sizes with a power-law we find the spin- 
stiffnes exponent A s = 0.056 ±0.11. The large error bar 
is due to poor statistics in the case of our largest size 
for which the sample average could only be taken over 
128 configurations of the bonds. The error in the deter- 
mination of the ground-state energy, estimated from a 
comparison of the results of searches conducted starting 
from different parent states, is much lower. It is worth 
mentioning that if we make the fit omitting the last point 
the result is A s = 0.052 ± 0.03. 

The smallness of |A| suggests that the LCD of the (± J) 
XY model is close to three. The data statistically favor 
A s > 0. If this is the case the spin-glass transition tem- 
perature should be finite with Tsg *C Tcg & s implied by 
our finding that A c S> |A S |. The numerical results are 
also compatible with a second and far less exciting pos- 
sibility, a zero-temperature transition with an unusually 
large correlation- length exponent. 

There exist no theoretical objections against a spin- 
glass transition below Tcg since the proofs 0,|l|l of the 
absence of SG ordering in the three-dimensional model 
at finite T fail if reflection-symmetry is broken ||[l2],[l3| 
as is the case when chiral-glass order is present. The 
possibility of a finite-temperature SG transition discussed 



here is thus restricted to the case of the XY model for 
which one can convincingly argue that Tcg 7^ 0. The case 
of the Heisenberg model seems quite different in that the 
role of the chiral variables is not as clear as in the planar 
model and there is no evidence of the existence of chiral 
order at finite temperatures. Therefore we still expect 
isotropic three-dimensional spin-glass models to have an 
ordered phase only at zero temperature |8p|,p^| . 
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